Hidden and antiferromagnetic order as a rank-5 superspin in URu2Si2 
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We propose a candidate for the hidden order in URu2Si2: a rank-5 E type spin density wave between Uranium 
5/ crystal field doublets F, ' and F^~ , breaking time reversal and lattice tetragonal symmetry in a manner con- 
sistent with recent torque measurements [R. Okazaki et al, Science 331, 439 (201 1)]. We argue that coupling of 
this order parameter to magnetic probes can be hidden by crystal field effects, while still having significant ef- 
fects on transport, thermodynamics and magnetic susceptibilities. In a simple tight-binding model for the heavy 
quasiparticles, we show the connection between the hidden order and antiferromagnetic phases arises since they 
form different components of this single rank-5 pseudo-spin vector Using a phenomenological theory, we show 
the experimental pressure-temperature phase diagram can be qualitatively reproduced by tuning terms which 
break pseudo-spin rotational symmetry. As a test of our proposal, we predict the presence of small magnetic 
moments in the basal plane oriented in the [110] direction ordered at the wave-vector (0, 0, 1). 



I. INTRODUCTION 

The nature of the low temperature ordered phase found in 
the heavy fermion compound URu2Si2 has defied explanation 
for over 25 years. While the transition into this hidden order 
(HO) phase appears quite conventional based on the effe cts on 
thermodynamic' ^, magneticP'^and transport phenomena^SMl, 
the order parameter itself remains elusive. This stands in con- 
trast to the significant entropy release^ across the transition, 
indicating a strong ordering, inconsistent with the weak anti- 
ferromagnetic moments that perplexed early studies. Heroic 
efforts have brought the full complement of experimental 
techniques to bear upon this problem, from local probes such 
as yuSR^ NMR** ^ and NQR'°, surface probe^pl scattering 
studies using elastic and inelastic neutrong^ISEl ^jj^j ^sonant 
X-ray^Mil] as well as explorations of the material through 
pressure, both hydrostatic'^ '^ '^ and uniaxial^" and high mag- 
netic fields^' . Key facets of the problem illuminated by these 
studies include a pressure induced antiferromagnetic phase 
(AF) and a superconducting phase that arises out of the HO 
phase. Similarities betw een the HO and AF phases, such as 
their Fermi surfaces^^ ^^ suggest a common underlying mech- 
anism. While great progess has been made, the central mys- 
tery of the identity of the HO still remains. 

Recently, an important clue to the nature of the HO has 
been seen in magnetic torque experiments^"*. Through a mea- 
surement of the off-diagonal magnetic susceptibility on small 
samples, it was found that the HO phase breaks the rotational 
symmetry of the crystal spontaneousl y in t he [110] direction. 
Lack of detectable lattice distortions^SE^ suggest this sym- 
metry breaking is purely an electronic phenomenon. This is 
inconsistent with many earlier theoretical proposals for the 
HO (sed^^ for a review) and thus has attracted a great deal 
of interest. New proposals to explain these observations in- 
clude, among other s, spin-nematic states^a dynamical sym- 



metry breaking --^yui staggered spin-orbit coupling ordei and 
hastatic ordei-'^H 

In this article, we propose a candidate for the hidden order 
in URu2Si2 as a rank-5 E type spin density wave between 5/ 



versal and the lattice point group symmetry D4;, in a manner 
consistent with the torque result. We argue that the expected 
coupling of this order parameter to magnetic probes, such as 
neutrons, can be effectively hidden by crystal field effects in 
the magnetic moment, while still contributing to second-order 
coiTelations such as susceptibilities. This would manifest as 
a small moment in the basal plane oriented along the [110] 
direction ordered at wave-vector (0,0, 1). In a simple tight- 
binding model for the itinerant heavy quasiparticles, we show 
the close relation between the HO and AF phases arises due to 
an approximate degeneracy between the E type (x and y com- 
ponents) and A2 type (z component) of the rank-5 pseudo-spin 
vector23 
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crystal field doublets F^ and F^ . This breaks both time re 



where Q = 2n/cz = (0, 0, 1 ) and we have defined gener- 
alized spin operators Sajs — \fa(^fi3 where a, (3 — 1,2 in- 
dicate the F*^' doublet. These doublets are related to the 
5/ 7o eigenstates as /i± - cos 6 f^i + sin 6* /^^i and /2± = 
- sin /j_ 5 + cos f^i where 9 is determined by the relative 

crystal field strengths^S The similarity of the Fermi surfaces 
arises as the breaking of the pseudo-spin symmetry is inop- 
erative in the k^t - ky plane, leading to similar cross-sections 
and thus oscillations. Using a phenomenological theory, we 
show the experimental phase diagram as a function of pres- 
sure and temperature can be qualitatively reproduced by tun- 
ing either the hopping which breaks pseudo-spin rotational 
symmetry or the effective coupling constants. As a conse- 
quence of this pseudo-spin symmetry breaking, the Goldstone 
mode connecting the HO and AF becomes gapped, leading to 
a resonance a finite frequency in the magnetic response. Spe- 
cific heat and magnetic susceptibilities are also presented, and 
experimental means to test our proposal are discussed. 

II. MODEL 

The heavy-fermion physics of URu2Si2 is signalled by a 
peak in the magnetic susceptibility^ indicating a coherence 
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FIG. 1. (Color online) (a) View of the Fermi surfaces in the ky^-k^ 
plane. The hole-like Fermi surface (red) is centered about fc- = Inic 
while the electron-like Fermi surface (blue) is centered about F. (b) 
View of the Fermi surfaces in the k^-k. plane. 



temperature ~ 60K, well above the hidden order transition 
temperature Tq- It is reasonable to assume that near the HO 
phase the relevant physics is encompassed by that of a (heavy) 
Fermi liquid. Recent work'^^'^"'" supports this, suggesting 
itinerant character for the U 5/ electrons, with coherent U 
dominant quasi-particles near the Fermi surface, can account 
for many of the properties of the paramagnetic state. 

Motivated by this, we consider single particle 5/ states in a 
crystal field of D4/, symmetry, using the notation of Ref.-''* for 
the irreducible representations of D4;,. The strong atomic spin- 
orbit coupling splits the 5/ levels into a low -lying j - Sjl 
sextet and high-lying j - 1 jl octet. Since the j - 1 jl states 
lie significantly higher in energy than the j - 5/2, we keep 
only the latter in our discussion. Under the potential of the 
crystal lattice this sextet spilts into three Kramers doublets, as 
r5/2 = IT-i © Fe, explicitly 

/i± = cos 61/^5 H-sin6l/:p3 

/2± = -sin0/^5 -Hcos6i/^3 

h± = f±\ 

We consider a tight-binding model for URu2Si2 using states 
with r^ ', F^ and F^ character, including contributions from 
nearest and next-nearest neighbor hoppings. The Ff, bands are 
involved in producing the incommensurate inelastic neutron 
pealP, at a nesting wave-vector connecting F7 and F(,, through 
the gapping of T-j Fermi surface in the ordered phases. For 
the sake of simplicity we focus only on Fermi surfaces of F7 
character as they form the superspin of the current study. If 
we include hoppings along the body diagonals, as well as in 
along the in-plane axes and diagonals symmetry restrict us to 
a Hamiltonian of the form 



Hi] - 2_^ [Alkfl^Jla,k + A2kfl^j2cr,k) + 
ka- 

X {Ckflj2.,k + Clflkf2-,k + h.c) + 
k 

X {Dkf'Ukh~,k - Dlflj2.,k + h.c) , (2) 

k 

where /„„., with a = 1,2, o" = + creates a state in the Fi, 
doublet with pseudo-spin cr. To obtain a Fermi surface that 



FIG. 2. (Color online) Graphical represenation of the relevant hop- 
ping amplitudes along the (a) body diagonal (b) in-plane axes and (c) 
in-plane diagonals shown in the unit cell of URu2Si2. 



matches our criteria above, we find it possible to take Ck = 0. 
The remaining dispersion functions take the form 

(akA i'^K\ (ak,\ 

Aak = 8fa COS — COS — I COS ^ I -H 

2t'g (cos (akx) + cos [akyj] + 

4tg cos (ak^) cos [akyj -1-1 + sgn(Q') , 
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where A12 is the eff'ective crystal field splitting between Fl, 

(2) 

and F^ levels, sgn(Q') = ±1 for a = 1,2 and c and a are 
the lattice constants. Physically, ta represents hopping along 
the a{x + y) + ^z and equivalent directions, t'^ along ax and 
ay and t'^ along a(x + y) and equivalent directions, as shown 
in Fig. [2] In the following we fix hoppings t\ - t2 - -0.3, 
t\ = -0.87, 4 = 0.0, t'( = 0.375, q = 0.25 with chemical 
potential ji - -0.5 and crystal field splitting A12 - 3.5. For 
concreteness we take |ri2| = 0.7 in the Fermi surface plots, but 
will allow it vary when discussing phenomenological mod- 
els. We note that Fermi surface does not depend strongly on 
the choice of fi2. Geometrically, one expects this f 12 hopping 
to be related to hybridization with the Ru d orbtials. As we 
have neglected several bands from our description, we treat 
the remainder of the system as a charge resevoir, keeping the 
chemical potential fixed and allowing the density vary. We 
note that for t\2 - this tight-binding model has global SU(2) 
pseudo-spin rotation symmetry within each doublet. For finite 
fi2 this is broken to a single extra U(l) symmetry between 
the F^ and F^ doublets, transforming /lo- — » e'"''^ f\a- and 

f2a ^ e--V2^- 

Using the above hopping parameters, the Fermi surface 
shown in Fig. [T] is obtained. These parameters have been 
chosen to respect the following experimental and theoreti- 
cal results. From Hall and magnetoresistivity measurements, 
one expects closed electron and hole Fermi surfaces of nearly 
equal size. Ab-initio calculations suggest that these electron 
and hole pockets are composed mainly of F^ and F^ . Fur- 
thermore, the similarity of the Fermi surfaces in the HO and 
AF states and nesting between the electron and hole Fermi 
surfaces at the AF ordering vector Q, strongly suggests that 



the HO and AF orderings occur at the same wave-vector. 
Under this assumption we then need the nesting to be im- 
perfect to match the small pockets observed in oscillation 
measurements'**^^. This is implemented in a fashion consis- 
tent with ab-initio results, with the electron pocket near F be- 
ing elongated in the [100] and [010] directions and the hole 
pocket at Z being elongated in the [1 10] and [1 10] directions. 
After folding along Q one then has four small pockets along 
[100] and [010] as shown in Fig|5] 



each site in the F7 basis 



ni. ORDER PARAMETER 



The key signals found in the torque measurements are 
Xxy ^ andy^j: - Xyy within the HO phase. Landau-Ginzburg 
arguments^ show that an order parameter consistent with 
these facts and their temperature dependence should transform 
as the two-dimensional E irreducible representation"*^ of the 
Dnh point group, with periodicity at the wave vector Q. The 
details of the torque oscillations single out a set of symmetry 
related directions in this two dimensional space, namely the 
four orientations (±1, ±1), so we will denote this type of or- 
der parameter as £(1, 1). In the picture outlined above, where 
the relevant degrees of freedom are the two F7 bands, there 
are only two local E type order parameters that mix /] and /2, 
taking advantage of the nesting at Q. These are the x and y 
components of the vectors / (5 12 - 021 ) or 5 12 + 021. 

The first E type order parameter is inconsistent with 
resonant X-ray scattering resultsl^l as it carries an electric 
quadrupole moment. The second breaks time-reversal and, 
as the in-plane components of the magnetic moment operator 
also transform as E, an induced magnetic moment is expected 
generically. This can be seen explicitly by restricting to the F7 
doublets and constructing the moment operator fl - L + 2S at 
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FIG. 3. (Color online) (a) The off-diagonal magnetic susceptibility 
Xxy in the £(1, 1) at a function of temperature. One can obtain a 
strong resemblance to the measured susceptibility iil^, shown in red, 
by fitting the overall amplitude of ^v^/_;^„ and the position of To- (b) 
The uniform magnetic susceptibility in the £(1,1) phase as a function 
of temperature. Both;^'^^- = Xm (red) and;^f,^ (blue) are shown. 



f^x^ 



6V5 

6V5 



[sin (26) [S [1 - S 22) + cos (26) [s [2 + S 21)] , 



f^y-^-[sm(20){s'{,-SQ + cos(20){sl,+Sl,)]. 
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[cos{20){s\,-Sl^)-sm(2d){s\,+Sl,)] 



6 r ^ 1 

+ 7 [5-11+^22] ■ 
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For 6-^+6 with 6 small, the terms proportional to 012 + 021 
in the x and y components are 0(6), while the z component 
remains 0(1). One has the scenario, that for small S the mag- 
netic coupling to this type of E order parameter is significantly 
reduced^^. For example the neutron scattering cross-section is 
0{6^). For these reasons we take the x and y components of 
Eq. [ijas the HO with the z component as the AF'^^. Note 
that even though the expectation is supressed, the presence of 
other terms ~ sin (26) in the moment operator allows terms 
of 0(1) to present in the in-plane susceptibilities. The bare 
susceptibilities in the £(1, 1) phase are shown in Fig. [3] The 
temperature dependence of the thexxy signal is strikingly sim- 
ilar to that found in the torque experiments. Deviations at the 
lowest temperatures are likely due to effects from the super- 
conducting phase. 



IV. PHENOMENOLOGICAL THEORY 

Due to the approximate nesting of the Fermi surfaces, the 
similarity to the AF state and the implications of the Landau- 
Ginzburg analysis of the torque experiments^'*, we consider 
only ordering at wave-vector Q. To explore this scenario, con- 
sider the phenomenological model, where the tight-binding 
model is supplemented by the coupling terms 



Hj, — — 



^■J]e<^'-{^n(i) + ^2di)), (4) 



4> 
'2V2 , 



where the momentum sum runs over the reduced Brillouin 
zone, folded along Q, as well as the quadratic terms 

where we have introduced the couplings gxy and gr. The 
£■(1, 1) phase cooresponds to \(f)x\ = |0^| = (p and 0. = 0. 
The phase where (px - (py - Q and (p^ t^ is magnetic, with- 
out any suppression factors of 5 in the magnetic coupling, and 
does not break the C4 symmetry present in the point group. 
We identify this with the pressure induced AF phase 

Explicitly, we first write the kinetic part of the Hamiltonian 
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where *}*]. = 



(/iV,* /i-,t fLk fl.k) and n = i(TyR&Dk + 
io-JroDi^. Then the phenomenological HamUtonian is given 
by 
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where the momentum sum now runs over the reduced Bril- 
louin zone. In this general form one cannot access the spec- 
trum analytically of Eq. |8] but since for k^ - the off-diagonal 
term D/^ is zero, one can find the spectrum in this plane, yield- 
ing branches fij and E^^„ where 
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This implies that the Fermi surfaces have the same cross 
sections at k, = 0, independent of the orientation of 
(p. This naturally explains the similarity in the oscillation 



measurements in the HO and AF phases and the remain- 
ing small pockets. The Fermi surface for various values of 
101 is shown in Fig. pi As the ordering is strengthened, the 
remaining Fermi surfaces form four small pockets along the 
X and y axes. This removal of large sections of the Fermi 
surface within the HO phase is qualitatively consistent with 
implications from transport and specific heat measurement^. 
The destruction of the pockets along the diagonal directions 
occurs due to more robust nesting in planes along the intersec- 
tion points of the electron and hole Fermi surfaces. A similar 
feature has been noted in recent ab-initio Fermi surfaces-"*. 

A key signature of the HO phase is the inelastic neutron 
scattering peak at Q. To examine this the imaginary part of 
the bare susceptibility at the ordering vector Q, ImxzziQ, co) is 
computed and shown in Fig. l4](a). There is a clear gap below 
a peak at finite frequency, reminiscent of the neutron results. 
The specific heat shown in Fig. l4](b) also shows a clear jump 
at To as expected, but the calculated value is too small to be 
realistically compared with the experimental results. The pre- 
cise size of jump in the model depends on the strength of order 
parameter, (p, tendency of nesting on the Fermi surface, inter- 
action effects and details of the some the other bands that have 
been left out of our model (for example from the gapping of 
the incommensurate modes). 

To study the transition for the HO to AF phases we look 
at the mean field phase diagram as a function of temperature. 
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FIG. 4. (Color online) (a) The specific heat as a function of temper- 
ature. Note the jump at To- (b) The imaginary part of the suscepti- 
bility in the zz channel at Q ai T = 0. The frequency w is scaled by 
the gap at perfect nesting A = \(p\/ V2. 



effective coupling constants g^y and g^ and t\2, the hopping 
which breaks pseudo-spin rotation symmetry. One can evalu- 
ate the free energy for this Hamiltonian numerically, including 
the quadratic parts of Eq. l5] and minimize to find the value of 
^. As the locking of the pseudo-spin orientation to the spa- 
tial orientation is controlled by the tn hopping, to obtain an 
£(1, 1) phase over a £(1,0) or £(0, 1) phase we choose the 
phase of ti2 to be n/A. This is supported by Slater-Koster type 
calculations of f 12, which also yield a phase of 7r/4. We find 
that for any finite value of f^ with equal couplings, g^y = gz, 
the AF state is favoured. This allows for a scenario where 
Oxy > 9z ^nd fi2 increases as pressure is raised from ambient, 
causing a transition from HO to AF as shown in Fig. l6](a). A 
more direct possibility is allowing gz/gxy to vary from g„j > gz 
at ambient pressure to gz > g^j at higher pressure, and fixing 
fi2, as shown in Fig. l6](b). In both these cases the HO to AF 
transition is first order, and the quaUtative topology is simi- 
lar to that seen as a function of pressure and temperature in 
URu2Si2. 



V. DISCUSSION 

An important ingredient in this proposal is the fine-tuning 
of the crystal field angle to be near enough to n/A to have the 
basal plane moments avoid detection thus far. Depending on 
the magnitude of the deviation 6, the order parameter should 
be observable in polarized neutron scattering. This leads to 
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FIG. 5 . (Color online) Fermi surface after imposition of finite (p = \ip\ 
in the k- = plane for several values of ^ (in the unfolded Brillouin 
zone). 
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FIG. 6. (Color online) (a) Mean field phase diagram as a function 
of the hopping parameter fi2, which breaks the HO-AF degeneracy, 
^t g-Jgxy = 0.8. (b) Phase diagram as a function of the coupling 
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5.0. 



a prediction of small moment in the basal plane in the HO 
phase, oriented along [110]. So far, neutron scattering has 
been focused on scattering at wavevector (1,0,0), favourable 
to detecting ordering oriented along the z axis. Furthermore, 
when (un-polarized) scattering is done at wave-vector (1, 0, 0) 
the small moment from the HO would be lost in the signal 
from the parasitic AF moment. To avoid a signal from the 
parasitic moment present in the HO phase, while still detect- 
ing a small moment along [1 10], we suggest a careful neutron 
scattering analysis with wave-vector (0, 0, 1) (parallel to Q) 
since the z component moment would not contribute due to 
the geometric factor in the scattering cross section. In par- 
ticular, when neutron is polarized parallel to g, all magnetic 
scattering is spin-flip and should come from the basal plane 
moment. Polarized scattering along (1,0,0) should also be ef- 
fective at uncovering the moment, so long as the sensitivity is 
sufficiently high. While the quantitative determination of the 
angle 9, and thus expected moment size to be seen in experi- 
ment, is beyond the scope of the current study recent LDA-hU 
calculations^^ suggest that such fine-tuning could be realized 
in URu2Si2. Effects of high magnetic fields could also lead to 
testable implications of this model. One expects that at large 



enough fields this will destroy the HO due to breaking of nest- 
ing between the F^ and F^"^ orbitals. The complex structure 
of the magnetic moment operator in Eq. l3] could lead non- 
trivial Fermi surface splittings, perhaps leading way to another 
type of ordering. The relation of this to the phases^' which ap- 
pear in the vicinity of ~ 37T warrants a detailed study, which 
we leave to future work. 

In summary, we have presented a simple effective tight- 
binding model for the relevant bands in URu2Si2. Motivated 
from this model, we propose that the HO is realized as a rank- 
5 inter-orbital spin-density wave oriented along the [110] di- 
rection. This gives the correct breaking of tetragonal sym- 
metry to account for the recent torque results^'*. A natural 
relation between the HO and AF phases is provided, connect- 
ing them through a pseudo-spin rotation as well as explaning 
the similarity of their Fermi surfaces. The Goldstone mode 
in the HO is gapped due to terms that break pseudo-spin rota- 
tional symmetry, manifesting as the finite frequency peak seen 
in inelastic neutron scattering experiments. The pressure- 
temperature phase diagram is also understood through this 
symmetry breaking. While the size of the moment in the basal 
plane in the HO can be suppressed by fine-tuning of the crys- 
tal field strengths, it should be detectable by polarized neutron 
scattering. 

Note added: After completion of this work, we have been 
made aware of a first principles LDA-hU studjGIl which has 
proposed an time-reversal breaking E-type order parameter as 
a strong candidate for the HO phase in URu2Si2. The small- 
ness of the moment in this work lends support to our current 
proposal. 
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